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Abstract 

We study electric field quench in A/” = 2 strongly coupled gauge theory, us¬ 
ing the AdS/CFT correspondence. To do so, we consider the aforementioned 
system which is subjected to a time-dependent electric field indicating an 
out of equilibrium system. Defining the equilibration time teq, at which the 
system relaxes to its final equilibrium state after injecting the energy, we 
hnd that the rescaled equilibriation time k~^teq decreases as the transition 
time k increases. Therefore, we expect that for sufficiently large transition 
time, k —)■ oo, the relaxation of the system to its final equilibrium can be 
an adiabatic process. On the other hand, we observe a universal behavior 
for the fast quenches, fc -C 1, meaning that the rescaled equilibration time 
does not depend on the hnal value of the time-dependent electric held. Our 
calculations generalized to systems in various dimensions also conhrm univer¬ 
salization process which seems to be a typical feature of all strongly coupled 
gauge theories that admit a gravitational dual. 
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1 Introduction and Result 

Quantum quench is defined as the process of (rapidly) altering a physical 
coupling of a quantum system. Ever since the quantum quench became 
accessible in the lab experiments [1], it has attracted a lot of attention to 
investigate its various properties. Another reason for its appeal is that, due 
to the time dependence of the coupling, the system is out of equilibrium and 
therefore understanding the evolution of the system is interesting enough. 
However, only with weak coupling assumption, analytic calculations have 
been done successfully [2]. 

We are interested in considering a strongly coupled system in which a 
time-dependent electric field (rapidly) changes from zero to a finite value, 
electric field quench. Since charged particles, such as quarks, exist in our 
system, applying a time-dependent electric field generates a current which 
non-trivially varies over time. The physical quantity we would like to study 
is equilibration time. It describes the time at which the out of equilibrium 
system, where its characteristic is a non-trivial time-dependent current, ap¬ 
proximately reaches its equilibrium specified by a constant current. In prin¬ 
ciple, this time depends on the final value of the electric field and on the 
transition time for which the electric filed reaches its final value. Because of 
the strongly coupled nature of the system, the perturbative methods are not 
applicable and thus the AdS/CFT correspondence (or more generally, the 
gauge-gravity duality) presents itself as a reasonable candidate. 
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According to the AdS/CFT correspondence type IIB string theory on 
AdSsxS^ background, resulting from near horizon limit of a stack of D3- 
branes, corresponds to A/" = 4 super Yang-Milles (SYM) theory with the 
gauge group SU{N) in 3 + 1 dimensions [3-7]. In the large t’Hooft coupling 
constant and large N limit, this correspondence reduces to a duality between 
a classical gravity (on asymptotically AdSsxS® backgrounds) and a strongly 
coupled SYM theory. In other words, the duality enables us to study various 
aspects of the SYM theories at strong couplings by using classical gravity. 
Applying this idea, quantum quenches in strongly coupled SYM theories 
with gravity duals have been extensively studied, for instance see [8-13] and 
universal behaviour for fast quenches has been reported in [14-17]. As well, a 
generalization of the AdS/CFT correspondence states that a strongly coupled 
SU(N) gauge theory in p + 1 dimensions is dual to a classical gravity on the 
near horizon geometry of Dp-branes [18] . 

On the gauge theory side, in order to get a non-zero current we need 
charged particles, similar to quarks in the quantum chormodynamics, which 
transform according to the fundamental representation of SU(N) gauge group. 
From the view point of the AdS/CFT correspondence, on the gravity side 
the fundamental degrees of freedom can be added to the system by intro¬ 
ducing Nf Dg-branes into the background, sourced by the Dp-branes, in the 
probe limit by which we mean Nf/N [19]. In fact, the probe limit guaranties 
that the back-reaction of the Dg-branes on the background is negligible. Dif¬ 
ferent aspects of these systems, i.e. Dp-Dg systems, have been studied, see 
e.g. [20,21]. The dictionary of the AdS/CFT correspondence states that the 
asymptotic behavior of the shape of the probe brane gives the mass of the 
fundamental matter and the expectation value of the corresponding operator 
in the gauge theory [6]. 

By choosing suitable configuration for the gauge fields living on the probe 
branes, one can consider a constant electric held on the brane and study 
its effect on the physical quantities such as the mass of the fundamental 
matter [22]. In an interesting paper [21], it was shown that this electric held 
gives rise to a constant current on the brane and this current increases in the 
presence of a'-correction [23]. Recently, the ehect of time-dependent electric 
helds, which clearly indicate out of equilibrium system, has been studied 
in [24] where it was shown that a time-dependent current is induced on the 
brane which hnally relaxes to its equilibrium value, in accordance with [21]. 

As we already mentioned, we consider a A/” = 2 strongly coupled SYM 
theory which is subjected to a time-dependent electric held. The gravity 
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dual is described by embedding a single (iVj = 1) probe D7-brane in the 
AdSsxS^ background when a time-dependent electric held is turned on the 
probe brane. We would like to study the effect of the transition time k and 
of the final value of the electric field Eq on the equilibration time t^q- Our 
results can be summarized as follows: 

• Our calculations indicate that a universal behavior emerges when /c < 1, 
i.e. fast electric field quenches. More precisely, by universality we mean 
the rescaled equilibration time k~^teq becomes independent of the final 
value of the electric field for small enough k. 

• For k > 1, i.e. slow electric field quench, the rescaled equilibration time 
decreases as one raises k at fixed final electric field value. In fact, it 
seems that for large enough values of k, the relaxation of the system to 
its final equilibrium is an adiabatic process. 

We then generalize the above configuration to the Dp-Dg systems and dis¬ 
cuss the effect of dimensions, in which the degrees of freedom live, on the 
equilibration time. We observe that both slow and fast electric field quenches 
qualitatively manifest a universal behavior for the systems we take into ac¬ 
count. 

In Appendix A, we briefly review and explain the way that a time- 
independent current can be found by using the gauge-gravity duality. In 
this paper we are working with the pulse function (2.17b), which doesn’t 
have a continues second derivative, in Appendix B we introduce a new pulse 
function with a continues second derivative and discuss its physical outcomes. 
Our observations indicate that the main results remain unaffected. 


2 Time-Dependent Current 

Adding the fundamental matter (quark) to the A/" = 4 SYM is equivalent to 
consider the probe D7-brane in the asymptotically AdS background [19]. It 
is well-known that switching on a non-dynamical electric field on the probe 
D7-brane leads to a constant current at non-zero temperature as well as zero 
temperature [21]. The effect of a non-dynamical external electric field on the 
system, when the mass of the quarks has been supposed to be zero, has been 
studied in [24] and it has been then generalized for massive quarks in [25]. 
In the following a generalization of above system to Dp-Dg systems will be 
considered. 
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To do so, we start with the following supergravity solution corresponding 
to near horizon limit of N coincident Dp-branes [18] 

+ dxp^) + + u^dQ.l_p), 

dQl_p = d9‘^ + sin^ Od^l + cos^ Odily-p-k^ 

where 

H{u) = {-y-P, = h"-^, C 01 ..P = H-\ (2.2) 

u 

The Dilaton held is represented by 0 and C'oi..p is a (p + l)-form coupled to 
the Dp-branes. The length scale R can be written in terms of the string scale 
Is = y~(x' and the string coupling constant gs = e'^°° 

K-" = (2.3) 

For the case p = 3, the above solution reduces to AdS^ x S^. Since there 
is no decoupling limit for p > 5, we consider only p < 4 cases. Apart form 
the case p = 3, the Dilaton held is not constant, indicating that the gauge 
theory coupling constant is running and consequently the gauge theory does 
not enjoy conformal symmetry [18]. 

According to gauge-gravity duality, a strongly coupled SU{N) SYM the¬ 
ory living on the (p -|- l)-dimensional worldvolume of N coincident Dp-brane 
corresponds to a supergravity on the background (2.1) in the large N limit. 
The supergravity approximation is only reliable when the curvature and the 
string coupling constant are constrained to be small. These two conditions, 
in terms of a dimensionless coupling constant Qeff, result in 

l<Pe//< AT^, 

where 

9eff = 

9ym = 27rgs{27rlsy~^. 

We refer the reader to [18] for more details. 

The low energy ehective action describing the dynamics of the probe Dg- 
branes in an arbitrary background is given by 


(2.4) 

(2.5) 


S — Sdbi + Scsi 

Sdbi = -Tq j det{gab + Bab + {2'Ka')Fab), 

Scs = rqj 
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Tension of the branes is ^ = 

Ramond Bab are 

9ab 
Bab 

M, N, ..{a,b,..) are nsed to describe the spacetime(worldvolnme) coordinates. 
The backgronnd metric Gmn is introdnced in (2.1) and Fab is the held 
strength of the gange held living on the Dg-branes. In the Chern-Simons 
part, P[..] is the pnll-back of the bnlk spacetime helds to the worldvolnme of 
the probe brane and denotes the (n+l)-form Ramond-Ramond potential 
conpled to Dg-branes. 

In order to compnte a time-dependent cnrrent on the probe branes, it is 
more convenient to rewrite (2.1) nsing the following change of coordinates 

p = Msin6*, 

; (2.8) 

a = ucosu, 

and we hnally have 

= H~^ {—dti^ + dXp) + (^dp^ + p^dVil -|- dcr^ -|- (T‘^dVt^_p_^. (2.9) 

Now let ns consider the following embedding for the Dp-Dg system 

t X\ .. X(i Xfi-\-\ .. Xp p (7 

Dp X X X X X XX (2.10) 

Dq X X X X XX , 

where g = fc -|- d -|- 1. Now we need to introdnce a snitable ansatz to stndy 
a time-dependent cnrrent on the brane. Frist we consider the massless case 
meaning that all the transverse coordinates mnst be zero. In the presence 
of an external time-dependent electric held on the brane, we expect a time- 
dependent cnrrent (J*). According to AdS/CFT dictionary, the cnrrent 
is dual to the U{1) gauge held on the brane. Therefore, in our ansatz, we 
include ^ 

Aa:{t,z) = — J E{t')dt'+ a{t, z), (2.11) 

where a new radial coordinate z = is applied. Near the bound¬ 

ary 2 ; —>■ 0, the gauge held asymptotically approaches [21] 

a{t,z) = aoit) + 2U{2-Ka'Y^‘^ ^ (2-12) 


: {2TTYgsll'^^ and induced metric Pab and Kalb- 

= GuNdaX^dbX^, 

= BMNdaX^dbX^. 
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( 2 . 13 ) 


where Af = rq(27r)^ [21] and it is then easy to see that 


{J^(t)) oc \im d‘^a(t, z). 


In the background we are interested in, the Kalb-Ramond held is zero 
and regarding ( 2 . 10 ) conhguration, the form held in the background does 
not couple to the Dg-brane. Therefore, the CS-part of the action and the 
second term under the square root in ( 2 . 6 ) do not contribute to the action 
describing the dynamics of the probe branes. Moreover, in order to have a 
stable background at zero temperature we consider supersymmetric systems 
where q = p + A,p + 2,p and k = 3,2, 1, respectively. For such conhgurations 
we have p — q + 2{k — 1) = 0. With the above assumptions the DBI action 
for the background (2.1) reduces to the following Lagrangian 


C oc - j{27ia'n^r[{dtA,y - (a,24.)2], 

= Z^\/w, 


where Sdbi = f and 


a = 


2(P - 7) 
p-5 ’ 


a 


q — 2p + 9 /P — 5 


( 2 . 14 ) 


( 2 . 15 ) 


A few values of a, [3 and 7 , depending on the choice of {p, q) system, has 
been listed in Table 1. 


Table 1: Parameters 


Case 

p 

(1 

a 

(9 

7 

1 

3 

7 

4 

-5 

1 

2 

3 

5 

4 

-4 

1 

3 

4 

6 

6 

-7 

1/64 

4 

2 

4 

10/3 

-3 

(3/2)10/3 


The equation of motion for derived from the Lagrangian (2.14) is 


dr. 




w 


-a 




w 


= 0 . 


( 2 . 16 ) 
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Our aim is to solve the above equation of motion for for arbitrary values 
of p and q using the numerical methods. To do so. we need to specify the 


time dependence of the external electric held. We choose two functions 



(2.17a) 


if t < 0, 

ii 0 < t < k, (2.17b) 

a t > k, 


and will discuss their similarities and differences in the next two sections. 
The hrst function shows that the electric held is zero at f = — oo and it 
reaches Eq ai t = +oo. However, (2.17b) is zero from —oo to zero. At 
f = 0 it starts changing to reach a maximum value at t = k and remains 
constant afterwards. This function provide a more realistic alternative for 
(2.17a), as its behavior gives a more convenient meaning to the concept of 
turning on an electric held and, after certain ’’hnite” amount of time, reaching 
a maximum value Eq. The function (2.17b) has been used to describe a 
time-depend meson melting in the AdS-Vaidya background in [26] and its 
generalization in the presence of an external constant magnetic held has been 
done in [27]. Notice that in both functions, fc is a measure to indicate how 
fast the functions can reach a maximum value. In other words, k represents 
the transition time from zero to hnite value. For instance, a small (large) 
value of k indicates a fast (slow) growth in the electric held. 


3 D3-D7 Case 


Properties of the (3, 7) system which is subjected to the time-dependent 
function (2.17a) were studied in [24]. According to (2.12) and (2.13), due to 
a time-dependent electric held, a time-dependent current is induced on the 
brane. At early times the system is highly out of equilibrium and the current 
reaches its equilibrium value asymptotically. The equilibration time can be 
dehned as the time at which the hnal value of the time-dependent current is 
equivalent to its equilibrated value. In order to obtain the equilibration time, 
as it was mentioned in the previous section, we will use the pulse function 
(2.17b). Since this function does not have a continuous second derivative, 
in Appendix B we introduce another pulse function which has a continuous 
second derivative and study its features. 


<r(i)> 


<r(i)> 






Figure 1: Numerical results for current as a function of time for different 
values of electric fields, transition times and pulse functions, 
top: Eq = 0.001 with k = 0.5 (left) or k = 1.5 (right), 
bottom: Eq = 0.005 with k = 0.5 (left) oi k = 1.5 (right). 


k'U. 

4 

3 

2 

1 

0 . 


eq 



Figure 2: The left panel shows the resacled equilibration time versus the 
transition time. The right panel indicates that (k~^teg)~^ goes to zero for 
smaller values of transition times, as we shall see for fast quenches. 
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As it was expected and is clearly seen from fignre 1, the amplitnde of 
the pnlse fnnctions (2.17a) and (2.17b), Eq, and the time of energy injection, 
k, control the large peak of the time-dependent cnrrent, for instance, at 
hxed k by increasing Eq the initial cnrrent peak rises. Similarly nnmber of 
oscillations increases when one considers higher Eq. On the contrary, as one 
increases /c, the valne of the initial large peak and the nnmber of oscillations 
decrease for a hxed valne of Eq. Moreover, onr investigations show that the 
independence of the hnal cnrrent from the choice of the pnlse fnnctions, is 
a nniversal featnre of the system, as it was expectable. On the other hand, 
different behavior of the system at early time, dne to different pnlse fnnctions, 
indicates that the system is sensitive to the manner of energy injection. As 
expected, different methods of injecting energy to the system wonld canse 
different reactions dnring this early time interval. 

We continne this section by stndying the eqnilibration time. In [24], it 
was discnssed that the appearance of an apparent horizon on the D7-brane 
signals the thermalization on the brane. In fact the earliest time that an 
observer sitting at the bonndary can see the apparent horizon formation in 
the bnlk is considered as the time-scale of the thermalization. More explicitly, 
a time-scale for the thermalization in terms of the location of the apparent 
horizon has been dehned which the anthors claim that it behaves nniversally; 
meaning that the time-scale of the thermalization is independent of k. 

Here, following [8,11,12], we dehne 


e{t) 


Ui,) 


( 3 . 1 ) 


where e{t) denotes the relative error for the cnrrent. is the hnal eqni- 
librated cnrrent that J^{t) relaxes to in the long time limit. Therefore the 
eqnilibration time can be considered as the time that {J^{t)) deviates less 
than 5% from its hnal eqnilibrinm valne, i.e. e(feq) < 0.05 and stays blew 
this limit thereafter. Notice that (J*^) is given by R{J^g) = {27ia'E)^ [21] 
for D3-D7 system (see also Appendix A for more details). We fonnd that 
the response of the system to the energy injection becomes diherent for two 
cases, k > 1 and k < 1. We will snmmarize the resnlts in the following. 


• Slow Quench 

In the case of A; > 1, onr nnmerical resnlts are shown in hgnre 2. It is 
clear that by increasing the valne of the transition time fc, the rescaled 


10 





(5J 



max((5J) 

0.06; 

0.05: 

0.04; 

0.03; . 

0 . 02 : 

0 . 01 ; 

0.00 ^ . . ,. ***!**• ., k 

0 5 10 15 20 25 


Figure 3: Left: 5J = \{J^{t)) — {Jeq)\ in terms of E{t) has been plotted for 
k = 1 (blue), /c = 3 (red), fc = 5 (green) and /c = 10 (brown) (top to bottom). 
Right: the maximum of versus k. 


equilibration time k~^teq decreases. In fact, when for example k ~ 
25, the equilibration time becomes comparable to k and therefore the 
system doesn’t have enough time to equilibrate during the transition. 
For smaller values of fc, the system needs more time to equilibrate after 
hnishing the injection of energy. 

A comparison between a time-dependent current and an equilibrated 
current, which in our unit is (Jg^) = has been plotted in figure 

3(left). One can see that 6J = \{J^(t)) — {Jeq)\, which is a measure 
of deviation from equilibrium, is noticeable for small k. However, for 
larger values of the transition time this deviation becomes smaller and 
hence one expects that 5 J —)■ 0 when the transition time goes to infinity. 
Therefore, It is expected that for sufficiently large value of k, i.e. /c —>■ 
oo, the equilibration time becomes much smaller than the transition 
time indicating that the degrees of freedom living on the brane have 
enough time to thermalize during the transition. As a result, in this 
limit the system passes through approximately equilibrated situations 
and consequently the above limit resembles an adiabatic limit. In the 
right picture the maximum value for 6J is depicted versus A:, which 
clearly shows the smooth flattening of the graphs in the left figure. 
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0.010 

0.005 


0.000 


I * 


0.00 0.02 0.04 0.06 0.08 0.10 


Figure 4: The inverse of the resacled equilibration time versus the transition 
time for Eq = 0.24 (blue), Eq = 0.27 (red) and Eq = 0.31 (green). 


• Fast Quench 

The behavior of the equilibration time in the region with fc < 1 is 
interesting enough to investigate separately. This behavior can be seen 
in the data shown in hgure 4. As seen in these hgures, when fc < 1, 
for three different values of the electric held, their equilibration time is 
not the same. By decreasing the transition time, one can see that the 
difference among equilibration times become smaller. For small enough 
values of /c, the hgure shows that the equilibration time for these three 
values of the electric held become the same or in other words it seems 
that the equilibration time is independent of the hnal values of the 
electric held in the case of fast quench and it is set by the transition 
time k, only. 

Gravitational collapse of a massive scalar held in a time-dependent back¬ 
ground which is asymptotically AdS spacetime has been studied in [12], In 
the held theory side, the mentioned conhguration corresponds to a thermal 
quench in a strongly coupled A/” = 4 SYM theory. Using two-point correla¬ 
tion functions of suitable operators and the evolution of the entanglement 
entropy of the above system, the equilibration time has been computed. In¬ 
terestingly, similar to our case, for slow quenches the rescaled equilibration 
time decreases with an increase in the transition time. In the case of fast 
quenches, a universal behavior is also reported in the literature, for instance 
see [12,14-17]. 

We also like to emphasize on a signihcant point here. In [12] one should 
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solve the equations of motion for the background metric and scalar held, 
simultaneously. In other words, the back-reaction of the scalar hied has 
been considered on the background and therefore the Einstein’s equations 
describe dynamics of the system. On the other hand, in our case we are 
working in the probe limit and we ignore the back-reaction of the probe 
brane on the background. Therefore, DBI action describes dynamics of the 
system in our case. Despite the fact that the actions and equations of motion 
are diherent, the slow and fast quenches of the scalar held and electric held 
behave similarly. One may conclude that the universality observed in these 
results happen is generic property of strongly coupled theories with gravity 
duals, similar to the universality of shear viscosity [ 6 ]. 


4 Generalization to Dj9-Dg 

Now one can question whether the dimensions, in which the fundamental 
matter or SYM theory degrees of freedom live, ahect the hnal physical results 
such as equilibration time and the universal behavior at small k. In order 
to hnd out the answer, we consider Dp-Dg systems where depending on 
the embeddings of each brane, the various degrees of freedom can live on 
different dimensions. According to AdS/CFT dictionary, in the case of D3- 
D7 system the SYM theory degrees of freedom, in the adjoint representation 
of the corresponding gauge group, as well as fundamental matter are living 
on the 3-1-1 dimensions. A generalization of AdS/CFT dictionary to Dp-Dg 
system has been extensively studied, for example see [20,21]. In short, for 
the supersymmetric cases, g = p + 4, g = p-|-2 and g = p, the SYM theory 
degrees of freedom live on p-|-l dimensions. However, the fundamental matter 
is conhned to a d-1-1-dimensional defect, see (2.10), where Dp- and Dg-branes 
coincide. In the following we consider three interesting cases. 

• D2-D4 System 

From (2.10), it is easy to see that the supersymmetric brane conhgu- 
ration for this case is given by 

t Xi X2 p ^2 O' D 3 

D2 X X X (4.1) 

D4 X X XX , 

where the SYM degrees of freedom live in 2-|-1 dimensional space-time. 
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However, in order to preserve supersymmetry at zero temperature the 
fundamental matter is confined to 1 + 1 dimensions. 

• D3-D5 System 

The brane conhguration (2.10) tells us that the supersymmetric brane 
conhguration is 

t Xi X 2 X 3 p fl 2 O' 

D3 X X X X (4.2) 

D5 XXX XX , 

where the SYM degrees of freedom live in 3 +1 dimensional space-time 
and the fundamental matter is conhned to 2 -|- 1 dimensions. 

• D4-D6 System 

In this case the supersymmetric conhguration of the branes is given by 



t 

Xi 

X 2 

CO 

p 



DA 

X 

X 

X 

X X 



(4.3) 

DQ 

X 

X 

X 

X 

X 

X 



It is clearly seen that although the SYM degrees of freedom live on 
4 -|- 1-dimensional space-time, the fundamental matter is conhned to 
3-1-1 dimensions. 

We repeat the calculations in the previous section for above cases. In fact 
the results we are interested inm are the behavior of the (inverse of) rescaled 
equiliration time in terms of the transition time. These results has been 
represented in hgure 5. They show that, at least for the cases we have 
considered here, the universal behavior for the fast quenches and adiabatic 
process for the slow quenches are qualitatively independent of the dimensions 
in which the SYM or fundamental degrees of freedom live. 
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Figure 5: Numerical results for the (inverse of) resacled equilibration time, 
(left, fast quench) right, slow quench, hgure, versus the transition time. For 
these hgures we consider e{teq) < 0.1. 

Top (D2-D4): left plotted for Eq = 0.05,0.03,0.02,0.01 (top to bottom) and 
right plotted for Eo = 0.01. 

Middle (D3-D5): left plotted for Eq = 0.05 (orange), 0.03 (green), 0.01 (blue) 
and right plotted for Eq = 0.01. Note that there is a lot of overlaps among 
different colors. 

Bottom (D4-D6): left plotted for Eq = 0.05,0.04,0.03,0.02,0.01 (top to 
bottom) and right plotted for Eq = 0.01. 
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A Constant Electric Field in Y)p-Y)q Scenario 

In order to find as {J^{t —)■ oo)) = one needs to find a solution 

for (2.16) with 


Foi = dtA^it.z) 

— Eo, 

(A. la) 

Eqz dzAxiti z'j 

= 0, 

(A. lb) 

Fiz = d^Ax{t,z) 


(A.lc) 


According to (A.la) and (A.lb), w has no time-dependency and therefore 
substituting the above equations in (2.16) yields 






w 


= 0 . 


(A.2) 


The equilirated current, up to a constant C, is dehned as [21] 


{Je,) 


y/w 


(A.3) 


Substituting expressions for dtA^ and d^A^ from (A.la) and (A. lb) in (2.14) 
leads to the following 


w = 


1 - 


1 7 ( 2 ™')^ 2 
-I- - RC ^ -^0 

7(^![AZ„-a-2/3/ lx\2^ 
C'^Ra ^ \^eq/ 


and by choosing R 


1 and 27ra' 


w = 


1 - 



one obtains 


^z°‘Eq 


(A.4) 


(A.5) 


To ensure non-negativity of tc, we identify the single (real) root of the nu¬ 
merator and the denominator which yields 


= (A.6) 

The above formula, of course up to a constant factor, expresses (Teg), the 
equilibrated current, in terms of Eq (which can be considered as the maximum 
of the imposed time-dependent electric held). For the four cases under study, 
we present the power dependency of (J^^) on Eq as follows 
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Table 2: {J^^) — Eq dependency 


D3-D7 

(Je,) « 

D3 - D5 

(Jeq) OC Eq 

D4: — D6 

«> oc 

D2-DA 

(J-> a: 


kt, 


-1 


■eq 


/C“T, 


eq 




Figure 6: Numerical results for the (inverse of) resacled equilibration time, 
(left) right hgure, versus the transition time (for pulse function in (B.l)). 
For these hgures we consider e{teq) <0.1. 

The left panel is plotted for Eq = 0.05, 0.04, 0.03, 0.02, 0.01 and the right 
panel is plotted for Eq = 0.01. 


B Different Pulse Function 


In this section we would like to investigate whether the behavior of the 
rescaled equilibration time depends on the pulse function. Thus we will 
discuss another pulse function which has continuous second derivative. Let 
us start with 


0 


if f < 0, 


E{t) = Eq { eaf - 15(f)" + 10(f)3 if0<t<k, 


1 


if t > k. 


(B.l) 


The electric held is zero on (— cxd, 0) and starts rising to reach a maximum 
value aX t = k. It is easy to repeat the computations in section 3 and hnd 
the numerical results for D3-D7 case. These results are shown in hgure 6. 
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